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Abstract 

In this letter, we study some evolution networks that grow with linear 
preferential attachment. Based upon some recent results on the quotient 
Gamma function, we give a rigorous proof of the asymptotic Mandelbrot 
law for the degree distribution pk oc {k + c) ' in certain conditions. We 
also analytically derive the best fitting values for the scaling exponent 7 
and the shifting coefficient c. 

Complex networks are now the joint focus of many branches of researchl^"^!. 
Particularly, the scale-free property of some networks attracts continuous inter- 
ests, due to their importance and pervasiveness '^^^l In short, this property 
means that the degree distribution of a network obeys a power law P{k) ex k~'^ , 
where k is the degree and P{k) is the corresponding probability density, and the 
scaling exponent 7 is a constant. A pioneering model that generates power-law 
degree distribution was presented by Barabasi and Albert (B A) . 

In recent studies, it was found that in some complex networks, e.g. trans- 
portation networks and social collaboration networks 1*^, the degree distribu- 
tion follows the so-called "shifted power law" P{k) (X {k + c)-^, where the 
shifting coefficient c is another constant. This property is also called "Mandel- 
brot law" [101. 

To understand the origins of such Mandelbrot law, Ren, Yang and Wang[iil 
proposed a interesting growing network that is generated with linear preferential 
attachment. In such networks, there exits a recursive dependence relationship 
between every two consecutive degrees 
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where where fc = 2, n, n is the number of nodes, m is a positive integer 
constant and /3 € [0, 1] is another constant. 

Defining a = — 1, 6 = , we can abbreviate Eq.(l) as 

Pk[k + b]=pk-i[k + a] (2) 

To derive the asymptotic of the degree distribution, Ren, Yang and Wang[i-'^1 
studied the following three kinds of approximations: 



1 



I) forward-difference approximation, assuming 



^ « Pik) - p{k - 1) = p{k) ^pik) = ^pik) 
dk k + a k + a 

we have an estimation of tlic power-law as 

p{k) oc {k + ay^'''"^ 
II) backward-difference approximation, assuming 

dp(k) ,, ^, fc + l + a ,,, ,,, a-b ,,. 

we have another estimation of the power-law as 

pik) cx (/c + 6+ 1)"^''""^ 
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III) Suppose we must have a Mandelbrot law p(fc) cx (fc -|- c)~'^. As a result, 
we have — 1) cx (fc — 1 -I- c)~'* . Substitute these two approximations in the 
logarithm type of Eq.(2), we have 
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Rewrite Eq.(7) as 
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and apply the second order Taylor expansion of ^ in Eq.(8), we have 

A b + a + l\-^''-''^ 
p{k) cx fc H — 
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All these three estimations indicates that the scaling exponent of the degree 
distribution should be —(6 — a). Simulation resultsl^^' show that Eq.(9) gives 
the best approximation accuracy of the empirical distributions. However, we 
still need a rigorous proof of this interesting finding. 

Indeed, further assuming Pi^) = 1; we have the following matrix equa- 

tion 
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Using Gaussian elimination algorithm, we can directly solve p{n) from Eq. (10) 



as 
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Based on the recursive relationship Eq.(2), for a given ?i, we have 
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where fc = l, ...,n— 1. 

It is well known that for Gamma function V{z), we have r(z + 1) = zT{z). 
So, we get 

r(j + i + fo) , , r(j + i + a) 
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where j — 1, n — 1. 
From Eq.(12), we have 
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where A = p{n) (jXj^i j+F) r(i+aj ^ constant. 

Eq.(14) indicates that p{k) has the same asymptotic behavior of r(fc+i+bj • 
Actually, the quotient of two Gamma functions is a difficult problem that re- 
ceived consistent attentionsl^^"^^! . There are numbers of approximation formu- 
las which are not accurate enough for the above applications. Fortunately, an 
important results had been obtained very recentlyl^^l as 

Lemma l^^^l Given two constants s and t, when x oo, we have 
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where Fk{t, s) are the polynomials of order n defined by 

Fo(t,s) = l 

1 y^(_i)fc+i '^'=+i(^) ~ Bk+ijs 
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where n > 1, -Bfc(i) is the Bernoulli polynomials (page 40 of [16]) for t. 

Based on Lemma 1, from Eq.(14), we can have an accurate expansion of 
the degree distribution as follows 
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i k + Thus, we reach the following 



As fc — >■ oo, we have 
conclusion rigorously. 

Theorem 1 The degree distribution follows an asymptotic Mandelbrot law 
Eq.(9) for some complex networks that grow with linear preferential attachment 
depicted by Eq.(2). 
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